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ON ; 

, Wc introduce functional methods to study the non-equihbrium dynamics of a quantum massless scalar field 

at finite temperature in a gravitational field. We calculate the close time path (CTP) effective action and, using 
its formal equivalence with the influence functional, derive the noise and dissipation kernels of the quantum open 
, system in terms of quantities in thermodynamical equilibrium. Using this fact, we formally prove the existence 

' of a fluctuation-dissipation relation (FDR) at all temperatures between the quantum fluctuations of the plasma in 

thermal equilibrium and the energy dissipated by the external gravitational fleld. What is new is the identiflcation 
of a stochastic source (noise) term arising from the quantum and thermal fluctuations in the plasma fleld, and 
the derivation of a Langcvin-type equation which describes the non-equilibrium dynamics of the gravitational fleld 
influenced by the plasma. The back reaction of the plasma on the gravitational field is embodied in the FDR. From 
the CTP effective action the contribution of the quantum scalar field to the thermal graviton polarization tensor 



O 



(N 

> , 

\^*^ can also be derived and it is shown to agree with other techniques, most notably, linear response theory (LRT). 

QiQ , We show the connection between the LRT, which is applicable for near-equilibrium conditions and the functional 

' methods used in this work which are useful for fully non-equilibrium conditions. 

(3 , PACS number(s); ll.lO.Wx, 05.40. +j, -04.62.+V 

00 ' 
ON 

^ , The behavior of a relativistic quantum field at finite temperature in a gravitational field has been studied before 

' by Gross, Perry and Yafi^e Gribosky, Donoghue and Holstein Q|, Rebhan and co-workers ||], de Almeida, Brandt, 
Frenkel and Taylor |^,|| for scalar and Abelian gauge fields. The thermal graviton polarization tensor and the 

■ effective action have been calculated and applied to the study of stability of hot flat /curved spaces, and dynamics 
of cosmological perturbations. Using linear response theory (LRT) [^-12|, Jackiw and Nair |^ derived the response 



I. INTRODUCTION 



function at high temperatures for a non- Abelian gauge field in a thermal plasma of charged particles. 

To describe screening effects and stability of thermal quantum gravity, one needs only the real part of the polarization 
tensor; but for damping effects, the imaginary part is essential. The gravitational polarization tensor obtained from 
the thermal graviton self-energy represents only a part (the thermal correction to the vacuum polarization) of the 
finite temperature quantum stress tensor, which under more general conditions (e. g. dynamical background) should 
contain also contributions from particle creation (arising from vacuum fiuctuations at zero and finite temperature). 
These processes engender dissipation in the dynamics of the gravitational field and their fluctuations appear as noise 
in the plasma. These latter two effects have not been sufficiently appreciated in most previous works in thermal field 
theory. 



* Present address 
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In this work we use open system concepts and functional methods a la Schwinger-Keldysh ||l^,|T^ and Feynman- 
Vernon |]l5| . This framework encompasses LRT which we illustrate in a calculation of the quantum corrections of 
the scalar field to the thermal graviton polarization tensor with results identical to that obtained before |3|Q| . More 
importantly, we derive the noise and dissipation kernels, prove that they satisfy a fluctuation-dissipation relation 
(FDR) P, p^ , p^ at all temperatures, and derive a stochastic semiclassical equation for the non-equilibrium dynamics 
of the gravitational field under the influence of the thermal plasma including all three aspects mentioned above. These 
cannot be obtained easily by LRT and other thermal field theory methods. 

Here the thermal plasma will be described by a massless scalar quantum field at finite temperature and the grav- 
itational field by a linear perturbation from Minkowski spacetime. Since the plasma is considered to be in thermal 
equilibrium it will be characterized by a time-like normalized vector field it^, representing the four- velocity of the 
plasma, and the temperature of the fluid in its own rest frame (3^^. Taking into account the four- velocity u'* of the 
plasma, a manifestly Lorentz-covariant approach to thermal field theory may be used However, in order to sim- 
plify the involved tensorial structure we work in the co-moving coordinate system of the plasma where = (I, 0, 0, 0). 

By making conformal transformations on the field and the spacetime, our results may be easily generalized to the 
case of a conformally coupled quantum scalar field at finite temperature in a spatially flat Friedmann-Robertson- Walker 
universe 0. Indeed we have earlier used the functional method and the Brownian motion paradigm ||l^ , p^ , |l8| , p^ 
to study similar problems in semiclassical gravity [po| . We found that quantum noise arising from fluctuations in 
the particle creation would constitute a stochastic source hitherto undiscovered [pf| (whose effect can overdominate 
the expectation value of the energy momentum tensor in the semiclassical Einstein equation j22'24|) in a new form 
of Einstein-Langevin equation p5|-|3^]. We also came to the understanding that back reaction of vacuum quantum 
field processes (such as particle creation) on the dynamics of the early universe near the Planck time is summarily a 
manifestation of a FDR in semiclassical gravity. 

The close time path (CTP) functional formalism of Schwinger [^3| is based on a generating functional defined in 
such a way that expectation values (rather than matrix elements) of time and anti-time ordered field products are 
directly computed. Keldysh developed a diagramatic technique approach useful for non-equilibrium processes 
in quantum field theory. This method has been further developed in [p^-^6[ (For some applications in field theory 
and cosmology see p7|- |4]] , p5| , ^ -p6[ ) . The advantage of this formalism is that effective dynamical equations for the 
expectation value of the quantum field may be derived from the effective action. Since these effective equations are 
real and causal a proper initial value problem may be formulated. The generalization of this formalism to curved 
spacetimes was given by DeWitt and Jordan, Calzetta and Hu p7| , ^ , ^ . 

In Sec. H, we describe our model and the derivation of the thermal CTP effective action. From this effective action 
we identify the dissipation and the noise kernels representing the linear response and the quantum fluctuations of the 
plasma, respectively, and obtain a stochastic semiclassical equation for the gravitational perturbations. This is done 
by using the underlying ideas of the influence functional formalism of Feynman and Vernon Sec. m is devoted to 
showing the existence of the FDR. We formally show that the noise and dissipation kernels extracted from the CTP 
effective action in general satisfy a formal relation at all temperatures. By explicitly computing the corresponding 
kernels we show that this relation is satisfied at high and zero temperature. In Sec. |^ we describe how our functional 
approach reproduces results from LRT. Finally, in Sec. 0, we give a brief conclusion. To make this presentation 
self-contained, in Appendix ^ we give a brief yet complete description of thermal field theory via the CTP functional 
formalism. 



II. CTP EFFECTIVE ACTION AT FINITE TEMPERATURE 
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A. The model 



In this section, we derive the CTP effective action for a thermal plasma in a gravitational field. To describe the 
plasma we consider a free massless scalar field 4> arbitrarily coupled to a gravitational field g^i, with classical action 



(1) 



where R{x) is the scalar curvature and the arbitrary parameter ^{n) defines the type of coupling between the scalar field 
and the gravitational field. If ^(n) = the quantum field is said to be minimally coupled, if ^(n) = (n — 2)/[4(n — 1)], 
where n is the spacetime dimensions, the field is said to be conformally coupled. To describe the gravitational field 
we consider a small deviation from flat spacetime 



(2) 



with signature (—,+,•••,+) for the Minkowski metric. Using this metric and neglecting the surface terms that appear 
in an integration by parts, the action for the scalar field may be written perturbatively as 

1 



1 hni/] — - 



d"x (j) 



(3) 



where the first and second order perturbative operators V^-^^ and T^'^^ are given by 

= ^{[d^h''''{x)]d, + ht'''{x)d^d, + an)R^'^\x)}, 
^2) = l\d^h'^-{x)]d^ + h^^x)d^d, - an){R^^^ix) + hix)R^'^ix)) 



(4) 



In the above expressions, i?^*^-* is the /c-order term in the pertubation h^i,{x) of the scalar curvature and the definitions 
hfj^^ and h^i, denote a linear and a quadratic combination of the perturbation, respectively, 

1 
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h 



For the gravitational field we take the following action 



(5) 
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d''x y^R{x) 



3R^,ap{x)R^''"''{x) - ( 1 - 360(e(n) - ] R{x)R{x) 



(6) 



4(n - 4) 

The first term is the classical Einstein-Hilbert action and the second divergent term in four dimensions is the coun- 
terterm used in order to renormalize the effective action. In this action £p = WwG, a = (28807r^)~"'^ and jl is an 
arbitrary mass scale. It is noteworthy that the counterterms are independent of the temperature because the thermal 
contribution to the effective action is finite and docs not include additional divergencies. 



B. CTP effective action 



In Appendix |a| we described how to compute the CTP effective action at finite temperature for a quantum scalar 
free field theory. In fact, we need to generalize the CTP formalism for two fields in the semiclassical approximation. 



This generalization, which is justified in Appendix AS , gives the following formal effective action (replace (j) by /i^j/ 



and V' by in equation (A24)) 
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r^TpihU = Sf^[h+,] - Sf^[h-^] - ^Tr{\nG^Jh^,]}, (7) 

where ± denote the forward and backward time path of the CTP formahsm and j, [h^^] is the complete 2x2 matrix 
propagator (a and b take ± values) with thermal boundary conditions for the differential operator □ + V^^'^ + V"*^^^ + • • •. 
In fact, the actual form of G^j, cannot be explicitly given. However, it is easy to obtain a perturbative expansion in 
terms of V^'^"*, the fc-order matrix version of the complete differential operator defined by = ^V^^ and = 0, 
and G^jj, the thermal matrix propagator for a massless scalar field in flat spacetime (A15-A16). To second order G^j, 
reads, 

_ T^(l)^/3 T/(2)^/3 , t a(1) t^(1) , /o^ 
'^ab - ^ab~'-'ac^cd ^ db ~ ^ ac^ cd db + ^ ac^ cd '^de^ef fb ^ V»j 

Expanding the logarithm and dropping one term independent of the perturbation h^^{x), the CTP effective action 
may be perturbatively written as. 



CTPl"'fJ.i^\ ^ ["-fii/i '-'g ["-fUJi 



.lTr[vi''^Gl^vi'^Gi^ + V^'^G^_V['''G^__ - 2vi'^ G^M''^ Gt^]. (9) 

The next task is to identify the noise and dissipation kernels from the non-local terms of the above thermal CTP 
effective action. The last trace is the only term that is responsible for dissipation because all others are constructed 
with symmetric kernels. In fact, this non-local term have information of the noise kernel also, as we will see in the 
next section when we show the existence of a FDR. In order to identify the role of each kernel we have to compute all 
the traces in the formal expression for the effective action. Some of these traces have divergencies that are canceled 
using the counterterms introduced in the classical gravitational action after dimensional regularization. In general, 
the non-local pieces are of the form Tr\Va^'^ G^^^vl^^'^ G^^]. In terms of the Fourier transformed thermal propagators 
^ab(^) these traces can be written as, 

d-xd-x' hl,{x)hip{x') J ^^e^H^~^')GL,{k + q)Giiq)T^-'^''{q,k), (10) 

where the tensor T^'''"^(g, k) is defined in Appendix In particular, the last trace of may be split in two different 
kernels N^"'-'°''^{x - x') and I)'"'-°''^{x - x'), 

'-Tr[vl^'^Gl_vi^^G^_+] = - f d'^xd^x' h+^{x)h-^{x')[r>^'''''f^{x~x')+iN^"''"^{x-x')]. (11) 



Recalling Eq. ( A16 ) one can express the Fourier transforms of these kernels, respectively, as 



N'^^'"''(fc) = TT^ / {9{k^ + q°)e{~q°)^e{-k°-q'')e{q°) + np{\q°\) + np{\k'' + q''\) 



-2n^(|g°|)n^(|fc° + q°\)} 5{qM{k + <z)']T'^'^'"'='(g, fc), (12) 



t'^'^'^^k) = -zn^l ^ {0(fc° + g°)0(-q°)-0(-fc°-g°)%'')+S9(fc° + g>^(|q°|) 

-sg(g°)nM|fc° -I- g°|)} S{q')S[{k + q)2]T'-'"/'(q, fc). (13) 

Using the property T^'''"''((j, k) = T^'''"^(-(7, -k), it is easy to see that {x- x') is symmetric and B'"''"'^ (x- x') 

antisymmetric in their arguments; that is, W^°''^{x) = N^'''"'^(-a;) and D^"'''^'^{x) = -B^"'''^'^ {-x). To properly 
identify these kernels with the noise and the dissipation of the system, respectively, we have to write the renormalized 
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CTP effective action at finite temperature (g) in an influence functional form jlj]. The imaginary part of the CTP 
effective action written in this form can be identified with the noise kernel and the antisymmetric piece of the real 
part with the dissipation kernel. This identification will be finally justified in two ways. In Sec. [II we will see that 
these kernels satisfy a thermal FDR and in Sec. EV we will see from the point of view of the LRT that D'"'-"'^{x) 
corresponds to the dissipation of the gravitational field or equivalently to the response of the thermal plasma to the 
gravitational perturbation, and N'^'''"''(a;) to the random fluctuations of the plasma. 

If we denote the difference and the sum of the perturbations h^^, defined along each branch C± of the complex time 
path of integration C by [/i^jy] = h~^i, — and {/i/iiy} = /i^^ + h~j^^ respectively, the influence functional form of the 
thermal CTP effective action may be written to second order in /i^^ as, 

1 



CTPl'^fiiy 



^^x d^x' [h^,]ix)L'(^f''{x - x'){h^p}ix') 



2£l 

+ ^ / d'^'x [h^,^]{x)T, 



d^x d^x' [h^,]{x)W^-'''l'{x-x'){h^p}{x') 
\ J d^xd^x' [h^,]{x)^^'''^''P{x-x'){Kp}{x') 



+ ^ / d^xd^x' [h^Mx)W'-^^{x-x')[Kp]{x') 



(14) 



The first line is the Einstein-Hilbert action to second order in the perturbation h^^{x). L'^^^"^{x) is a symmetric 
kernel (i.e. L''^^^°'^{x) — L'^^-j'^^ {—x)) and its Fourier transform is given by 



L(o) [k) 



A:2Tr"^(q, k) + 2fc"Tr"''(g, k) + T^'^iq, k) - 21:''^^''"' {q, k) 



(15) 



The fourteen elements of the tensor basis T^'''"'^(g, k) {i — 1, - ■ ■ , 14) are defined in Appendix In the second line 
Ty^-^ has the form of a perfect fluid stress-energy tensor 



30/?4 
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(16) 



where u^^ is the four-velocity of the plasma and the factor is the familiar thermal energy density for massless 
scalar particles at temperature /3~^. In the third line, the Fourier transform of the symmetric kernel II^'''"^(a;) can 
be expressed as 



^"^^"''(fc) 



tk^ ri 



ln^Q'-'"^(fc) + ^Q'^'''"'^(fc) 



1 ^^y.cit(3 1 

3" 



+ W {-Ti''^"''^^' ^) - 2Tr"^(«, k) + TT-\u, k) + 2Tr"''(u, fc)} 



+ ^ {A:2Tr"''(u, k) - 2fc2Tr"^(zi, k) - Tr"''(u, k) + l^^^^^^u, fc)} 



d'^q 



(fc + qf 



+ 8{{k + qf\n(i{\k° ^ q°\)V 



|T^"'"'3((7,fc), (17) 



where /i is a simple redefinition of the renormalization parameter /i given by /i = /iexp(Y| -I- ^ ln47r — ^7), and the 
''^(fc) and ( 

Q^"'"''(yfc) 



tensors Q'^'''"'^(fc) and Q'^'^'"^ (k) are defined, respectively, by 



T 



.[1 _ 360(C - 1)'] |Tr"''(g, k) + ^Tl!^'"^iq, k) - ^Ttr^iq. fc)| , 



(18) 
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= [1 + 576(e - -r 60(e - ^)(1 - 36C')] |Tr"''(g, fc) + ^T^^^'-'^lg, fc) - ^T^s'^^iq, fc) j ■ (19) 

In the above and subsequent equations, we denote the couphng parameter in four dimensions ^(4) by ^ and conse- 
quently ^' means d£^{n)/dn evaluated at n — 4. h'^ ' ^ [k) is the complete contribution of a free massless quantum 
scalar field to the thermal graviton polarization tensor ||^,^ and it is responsible for the instabilities found in flat 
spacetime at finite temperature [|l|,^,D. Eq. reflects the fact that the kernel '^(/c) has thermal as well as 
non-thermal contributions. Note that it reduces to the first term in the zero temperature limit (/3 — > oo) 

r""''(fc) ^ |^ln^Q--"''(fc) -f iQ''-"''(fc)} . (20) 

This asymptotic contribution to the CTP effective action includes both pure vacuum quantum fluctuations p8| and 
non-conformal fluctuations pO.O. On the other hand, the leading term (/3~"^) at high temperature may be written 



as 



H'""'(fc) - ^i:H.(.)Tr-''(u,if), (21) 



=1 



where we have introduced the dimensionless external momentum K^^ = = {r,k). The Hi(r) coefficients were 

first given in ||^ and generalized to the next-to-leading order (/3^^) in [Q. They have been reproduced with our sign 
conventions in Appendix |^ for completeness. It is important to note that the addition of the contribution of other 
kinds of matter fields, even graviton contributions, does not change the tensor structure of this leading contribution 
and only the overall factors are different j^]. However, this is not true for the subleading contributions as was shown 
in|]. 

Finally, the noise kernel W^^f^ix) and the dissipation kernel {x) are new and may be directly read off from 

( p^ ) and (|l3|), respectively. As expected the noise kernel appears in the imaginary part of the infiuence functional form 
of the CTP effective action and the dissipation kernel is the only antisymmetric kernel of the real part. Identification 
with the influence functional [ p^ provides a physical interpretation: N'^'''"^ {x) represents the random fluctuations of 
the plasma and T>^'^'°'^{x) the dissipation of energy of the gravitational field. 



C. Einstein-Langevin equation 



In this section we show how a semiclassical Einstein-Langevin equation can be derived from the previous thermal 
CTP effective action. This equation depicts the stochastic evolution of the gravitational field under the influence of 
the fluctuations of the thermal plasma. 

We first introduce the influence functional T = eicpiiSip) which is defined from the influence action 5*/^ of Feynman 
and Vernon |l^. In fact, we can write this influence functional in terms of the the CTP effective action because in 
the semiclassical limit this effective action is completely equivalent to Sj f [|i^,^,^ , 

^ = cxp* {Reivl'cTplhU} + ^Jd^^ d^x' [V](^)N^''"''(^ " x')[h^p]ix')^ , (22) 

where Re{ } denotes taking the real part. Following |l5y25| -|2^, we can interpret the real part of the influence 
fimctional as the characteristic functional of a non-dynamical stochastic variable j^'^{x), 

$([V]) - exp (^-i j d^x d^x' [h^,]{x)W''''P{x - x')[Kp]{x')^ . (23) 

This classical stochastic field represents probabilistically the quantum fluctuations of the matter field and is responsible 
for the dissipation of the gravitational field. By definition, the above characteristic functional is the functional Fourier 
transform of the probability distribution functional V[j^'^] with respect to j^". 
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([VD = /p/"' T'bn e'/"^'" '''""'^"^^'""^"^ (24) 



Using (23) one can easily see that the probabiUty distribution functional is related with the noise kernel by the formal 
expression, 

exp (-i/d^x d^x' j^Ax^'^^'^'^ix - x')]-^j„0{x')) 

nfl = ^ 7 T- (25) 

JVJ^^- cxp [-Ud^x rfV j^.(a;)[N^''^"''(x-a;')]-ija/3(x')) 

For an arbitrary functional of the stochastic field ^[j^''], the average value with respect to the previous probability 
distribution functional is defined as the functional integral {A[j^'^])j = J 'D\j^'^] V[j^'^]A[j^'^]. In terms of this 
stochastic average the influence functional can be written as !F ^ (exp (iF^yp[/iJ^]))j , where F^j,p[/iJ^] is the 
modified effective action 

^CTp[hU ^ Re{T^cTp[hU} + J d'x [h^.]ix)r'' (x). (26) 

Clearly, because of the quadratic definition of the characteristic functional ( p3| ) and its relation with the probability 
distribution functional (|2^), the field j'^'^{x) is a zero mean Gaussian stochastic variable. This means that its two-point 
correlation function, which is given in terms of the noise kernel by 

{nx)r''{x')), = w-^^^ix-x'), (27) 

completely characterizes the stochastic process. The Einstein-Langevin equation follows from taking the functional 
derivative of the stochastic effective action ( ^6| ) with respect to [/i^^](a:;) and imposing [hf^,j]{x) — |^ , ^ , ^ . In our 
case, this leads to 

^ J d^x' h1^f^ix-x')h^f3ix') + ^T^;^+ Jd^x' (R^''^"''ix-x')^B'^-'-''''ix~x'))h^p{x')+rix) = 0. (28) 

To obtain a simpler and clearer expression we can rewrite this stochastic equation for the gravitational perturbation 
in the harmonic gauge h'^'^'^^ = 0, 

ah^^ix) + fp I^Tl"^] + 2Pp,^^p J d^x' (H^'''"''(a: - x') - D^'''-'^\x - x')) hP^x') + 2f''{x)^ = 0, (29) 

where we have used the definition for h^'^{x) written in (^ and the tensor Ppa.ap is given by 

Ppa,afi = 2 iVpaVcrp + VpPVcra - VpcrVap) ■ (30) 

Note that this differential stochastic equation includes a non-local term responsible for the dissipation of the gravita- 
tional field and a noise source term which accounts for the fluctuations of the thermal plasma. They are connected 
by a FDR as described in the next section. Note also that this equation in combination with the correlation for the 
stochastic variable (^^ determine the two-point correlation for the stochastic metric fluctuations (Jif_iij{x)hai3{x')) j 
self-consistently. 



III. FLUCTUATION-DISSIPATION RELATION 



Now we want to see how the kernels found above to represent the noise and the dissipation of the system are 
functionally related. This FDR reflects the balance between the quantum fluctuations of the thermal plasma and 
the energy loss by the gravitational field. First, we explicitly see how this relation appears at zero temperature 
and compare with the results obtained in Refs. [ p8|j30| . Secondly, we discuss in some detail the relation in the high 
temperature limit. Finally, using the properties of the thermal propagators, we formally show that the FDR is satisfied 
in general. 
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A. Fluctuation-dissipation relation at zero temperature 

Since the thermal functions 't./3(|(7°|) vanish at zero temperature, the complete zero temperature contribution to the 
noise and dissipation kernels can be directly read from ([l^ ) and (|3|) 

^'T^{k) = ^'J^ ^il'Mk + q f] Mk" + qnOi-q") + - T'^^^'^^iq, fc), (31) 

f^iT^ik) = -^7r' J -0^ Siq^)6[ik + q)'] {^(F + - 9{-k" - T'^'^'^^iq, k). (32) 

The computation of these integrals is very simple because of the presence of the delta functions. Furthermore, since 
they are independent of the temperature the final tensor structure cannot depend on the four-velocity that defines 
the preferred co- moving coordinate system of the plasma. In fact, this simpler structure allow us to write all the 
integrals in terms of just only two integrals 

/ (1^ '^('^'^^[(^ + + ^ - = ^^"^'^%°!)t^^"^°^^ ' (33) 

one corresponding to the noise kernel (+ sign) and the other to the dissipation kernel (— sign). The final result for 
both kernels can be written in terms of the tensor (^^'^'"'^[k) defined in ( [T^ ) as 

&;:f'{k) = sg{k^) (^^) Oi^eW'^^'^f'ik). (35) 

It is clear from these expressions, that the FDR satisfied by the noise and dissipation kernel at zero temperature is 
given by 

Nr:r'(fc) = ~^ sgikni^',:r\k). (36) 

This zero temperature result includes in one expression the FDR found in p8| for conformally coupled fields and that 
found in for non-conformally coupled fields. We next discuss the pure thermal fluctuation contribution in the 
leading order at high temperature. 

B. Fluctuation-dissipation relation at high temperature 

Here, we describe in some detail the computation of the noise and dissipation kernels at high temperature in order 
to find the corresponding FDR. The computation of thermal integrals involving the thermal factors n/3(|g°|) is further 
complicated by the tensor structure of gravitation. 

At high temperature one considers that the external momentum is much less than the temperature; that is k°, \k\ <ti 
■ Here we assume that the temperature is high enough so that this condition is satisfied, but remains well below 
the Planck temperature <^ (ip^ ^ £p so that graviton loop corrections to the effective action may be neglected. 
On the other hand, the leading contribution can be obtained for internal momentum \q\ of the order of the temperature 
which means that the external momentum |fc| is small compared with the internal momentum i.e., |fc| ^ \q\ 
(see, for example ^^). Let us start with the leading contribution (/?^) for the noise kernel. The thermal contribution 
to the noise kernel n|^^) ^ (k) is given by the last three terms of Eq. ( p^ ) 
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Nf;f ''(fc) ^ n^J0^ Siq')S[ik + qf] {n^kl) + npHk" + g"!) + 2nMk1)r.M|fc° + 9°!)} T^'''"''(9, k). (37) 

Integrating the q° variable, and taking into account the properties of the tensor T^"^'"^ (q, k) the above expression 
simphfies to 

Nr;r'w = 1/(0 ^^[«Mo+«Mp°)+2n,(9>Mp°)] 

X [S{p° ~q°- k°) + 5{p° + g° + fc°)]T^'^^"'3(g, k) 
+{k ^ -fc}, (38) 



where (7° = |g| and p° — |fc + g|. On dimensional grounds it is clear that the leading contribution at high temperature 

-3 



to the tensor T^"'"^ {q, k) is given by its first term Ti^'^'"'^ {q, k) only. We are also going to assume the following 



approximations: 

p° ~ q'' + kq, 
p°-q°-k° ~ -k° + kq, 

p° + q° + k° ~ 2q°. (39) 
Further simplifications arise because in this limit the thermal factor n^(p°) may be approximated by 

n^(p°) ^ „^(g°) + fcq^!^. (40) 

For the particular combination of thermal functions in the noise kernel this yields 

np{q°) + np{p°) + 2np{q°)np{p°) ~ 2np{q°)[\ + np{q°)]. (41) 

As can be seen by a direct calculation the integrals with 5{p° + q° + k°) do not contribute in this approximation. 
Then, the evaluation of the leading contribution at high temperature for the noise kernel ( p8| ) is reduced to calculate 
a solid angle integral and standard thermal integrals, separately, 

NrJ^^W = TT^, I dnSiKQ)Ti;-'^^{Q,K) r dq° (g°)%(q°)[l + nMg°)]. (42) 

In order to make the solid angle integral dimensionless we have used the external momentum K defined previously 
and we have introduced a new lightlike vector = {l,q). Note that the theta function 6{k'^) appears as a condition 
imposed by the delta function in the solid angle integral. The explicit computation of these last two integrals in terms 
of the tensor basis T^'''"^(u, K) gives the final result to leading order in /3 

N^rW - 5$5^EN.MTr-'^("'^)' (43) 
where the 14 coefficients Ni(r) are given by 

8 

-5{K^ - 6r^)K'^, -5{3K^ - 4r^)K^r, iK^ - 24:K\^ + 8r'^,-{K^ - Ar^)K^, 5K^r] . (44) 



The determination of these coefficients involves the evaluation of the inverse of the cumbersome 14 x 14 matrix 
Mij = Ti ^jy^apiu, K)T'^'^'"^ {u, K) 1^,^. However, this can be easily performed with an algebraic computer program. 

The thermal contribution to the dissipation kernel i^^p) ^ {k) can be obtained similarly. From the thermal terms of 
Eq. (P) we have 
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Dl'jr^'Cfc) = -*7r2 J ^ S{q')d[ik + qf] {sg(fc° + - sg{q")npi\k'' + T^^^^^^iq, k). (45) 

This integral can be simplified to 

-5{p° + 9° + fc°)[n;3(g°) + n^(p°)]} T'^"^^"^(g, fc) 
-{fc^-A;}. (46) 

Again the terms with 5{p° + q° + k°) do not contribute to the final result. Using ( |40| ) to approximate the difference 
of thermal functions and following the above computation for the noise kernel, we can write the leading contribution 
in [3 to the dissipation kernel as 

Dr;r'w - -^^EN.MTr-'^i-^i^), (47) 

where Ni(r) are the coefficients given in (|4|). It is clear that these leading contributions for the noise and dissipation 
kernels satisfy a FDR at high temperature of the form 

" ' [w) ^^^^ 

It is important to note that the leading contribution for the noise kernel depends on whereas that for the 
dissipation kernel depends on /3~^ and this is a direct consequence of the FDR. In contrast to the previous zero limit 
temperature FDR, this result may be interpreted to be a consequence of pure thermal fluctuations. We next show 
that formally there is a FDR valid to any order that includes both quantum and thermal fluctuations. 



C. Fluctuation-dissipation relation at finite temperature 



The existence of a FDR at finite temperature between the noise and dissipation kernels we have identified in 
the previous sections is very easy to prove under general conditions. We begin by showing that the FDR naively 
appears at the level of propagators as a direct consequence of the KMS relation 1 11 1 (see Appendix 0) . Then, using 



a generalization of this KMS relation, we see how the FDR is also satisfied by our noise and dissipation kernels. 

To obtain the FDR at the level of propagators we need to introduce the Schwinger and the Hadamard propagators. 
These propagators are defined as the thermal average of the anticommutator G{x — x') = —i{[(t>{x)^ 4){x')])p and the 
commutator g''^^ {x — x') = — i({0(a;), (j){x')})p, respectively. The first represents the linear response of a relativistic 
system to an external perturbation and the second the random fiuctuations of the system itself ||ll|,|l2| . Since we can 
write the KMS condition satisfied by the propagators and G^^_|^ in Fourier space as 

&_4k) = e^^°G^^_(fc), (49) 

Schwin 

example, in terms of G^J__(fc) alone 



the Fourier transform of both the Schwinger G{k) and the Hadamard gL^^ (fc) propagators can be expressed, for 



G{k) ^ &^{k)~&i_{k) ^ (e^'^" ^l)Gl_{k), 
Gf{k) = G\{k) + Gl_{k) = (e/''" + l^Gi_ik). (50) 



The FDR satisfied by these propagators follows inmediately from the above equalities pT|jr 
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G';>{k) = coth(^YjGW- (51) 

Obviously, this relation can also be recovered if we write the explicit expressions for the Fourier transform of the 
propagators 

G{k) = -2TTi sg{k°)S{k^), 
G^p\k) ^ -2Tn coth(^^^^S{k^). (52) 

To use this last approach in our case could be a very difficult task because one needs to compute the integrals for the 
noise and dissipation kernels explicitly. On the other hand, if we follow the first technique wc only need to generalize 
the KMS condition of Eq. ( ^9| ) to product of two propagators. This generalization reads 

& + {k + q)&+_{q) = e^'=°G^_(fc + g)G^+(Q), (53) 

and can be used to deduce the following formal identity 

il^ (i,^r,\riP („\^ril^ {u^„\n0 — ( \ \r<P ii,^r,\n0 (^\_riP /u^„\nf^ 



G'\{k + q)G'l_{q) + Gl_{k + q)G^_+{q) = COth\^^ J [G>'_^{k + q)G'l^{q) - Gl_{k + q)G%{q)]. (54) 

Finally, one only has to write, from the trace of Eq. (^0|) and the definitions (12) and (p^), the noise and dissipation 
kernels in terms of the propagators G£::p, respectively, as 



N''''"''(fc) - -i/(04 [G%(fc + <z)G![_(<z) + G|^_(fc + g)G%(g)]T^— '^(<Z,fc), (55) 



D"""''(fc) - i/(04 [G!+(fc + g)G^_(9)-G^-(fc + 'z)G!+(g)]T''-"^(g,ft), (56) 
and use the formal equality ( p^ to prove that they are related by the thermal identity 

^i"''''\k) = zcoth(^^)D"'''"''(ft). (57) 
In coordinate space we have the analogous expression 

W^'^^ix) = J d-'x' KFD{x-x')B^'' "f^ix'), (58) 
where the fluctuation-dissipation kernel Kpoi^ — x') is given by the integral 

KM^ - :r') - * / ^''^■("""'^ coth ( . (59) 



(27r)4 ^ V 2 

The proof of this FDR at finite temperature is in some sense formal because we have assumed along the argument 
that the integrals are always well defined. Nevertheless, the exact results obtained for the zero and high temperature 
limits indicate that the noise and dissipation kernels are well defined distributions ||5C| ]. The asymptotic analysis has 
also been useful to determine the physical origin of the fluctuations. 

IV. LINEAR RESPONSE THEORY 

The purpose of this section is to understand the connection between the LRT [p[- plj5l]j5^ and the functional 
methods we have used here. In the spirit of LRT the gravitational field is considered as a weak external source which 
imparts disturbances to the plasma whose response is studied to linear order. 
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Let us first recall the main features of LRT. Consider a system described by the Hamiltonian operator Ho initially 
coupled to an external driving agent linearly, say Aa. Since we are only interested in how the system responds to 
the external agent, and not the details of the agent, we will ignore the Hamiltonian for the external perturbation but 
write the complete operator Hamiltonian of the sytem as 

H = Ho + AJ'', (60) 

where J" is the current operator associated with the external agent. If the system is in thermal equilibrium before 
the external source is applied the first order response of the system to this external force is given by the thermal 
expectation value of the commutator of the current operator over its thermal average 

{[j^ix) - {J^ix)),, r{x') - {r{x'))p])p. (61) 

In contrast, the intrinsic quantum fluctuations of the system are described by the thermal average of the anticom- 
mutator. In our case, the conserved current operator is given by the stress-energy tensor T^^{x) as derived from 
the classical action. The objective of this final section is to show that the response and fluctuation functions for the 
stress-energy tensor considered in the LRT are equivalent to our dissipation and noise kernels, respectively. 
First, we write the classical action for the matter field to linear order in the gravitational perturbations, 

1 
2 

with the stress-energy tensor given by 

Note that T'^'^{x) is conserved if the classical unperturbed equation of motion for (p is satisfied and it reduces to 
the stress-energy tensor for a scalar field in fiat spacetime if ^ — 0. Alternatively, we can write the Hamiltonian 
formulation of our problem. If we introduce the conjugate momentum variable of the matter field to first order in the 
perturbation 

dC -1^ OT^"" 
2 



Sm[(f>,h^.] - 77 / d^x [000 + Vrn> (62) 



H EE — ^ <^+-V-7^, (64) 



the Hamiltonian can be written as 

H ^ \jd'x [tf + (V0)2 - /i,,,T'^'^] . (65) 

Note that to first order cf) and H are interchangeable in the expression for the stress-energy tensor. 

Using the thermal version of the Wick theorem [^ , ^ , one can write, after some algebra, the equilibrium thermal 
average of the two-point function for the stress-energy tensor operator at different spacetime points in terms of 
products of thermal propagators, 

{T^'^{x)T-P{x'))p-{T^'^{x))p{T"^{x'))p^-2 J 1 ^ &^ik + q)Gl_{q)T^'^^-^ {q, k). (66) 

Finally, defining ApT^'^{x) = T^'^{x) — {T^'^{x))p and using the expressions for the the noise and dissipation kernels 
given in (^5|) and ( p6| ) respectively, we obtain 

{{ApT^'^{x),A0T''0ix')})p = 8 W'^'^'' ix-x'), (67) 

{[ApT'^^x),ApT''^{x')])p = SiB^-^-t'ix-x'). (68) 

From these formal identities we conclude that the functional method gives a description of the lowest order dynamics 
of a near-equilibrium system equivalent to that given traditionally by the LRT. 
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V. CONCLUSIONS 



In this paper we show how the functional methods can be used effectively to study the non-equilibrium dynamics of a 
thermal quantum field in an external gravitational field. The close time path (CTP) effective action and the influcence 
functional were used to derive the noise and dissipation kernels of the open system. We show the formal equivalence 
of this method with the traditional linear response theory (LRT) for lowest order perturbance of a near-equilibrium 
system, and how the response functions such as the contribution of the quantum scalar field to the thermal graviton 
polarization tensor can be derived. An important quantity not usually obtained in LRT approaches but of equal 
manifest importance in the CTP approach is the noise term arising from the quantum and thermal fluctuations in the 
plasma field. With this we derived a Langevin-type equation for the non-equilibrium dynamics of the gravitational 
field under the infiuence of the plasma. The back reaction of the plasma on the gravitational field is embodied in 
a fiuctuation-dissipation relation (FDR), which connects the quantum fluctuations of the thermal plasma and the 
energy dissipated by the external gravitational field. We formally prove the existence of such a relation for thermal 
plasmas at all temperatures. 

This powerful method perfected in recent years has been successfully applied to semiclassical gravity and cosmolog- 
ical problems by a number of authors. It should also prove useful for thermal field theory problems in QED and QCD, 
such as electroweak phase transition, quark-gluon plasma and disoriented chiral condensates in heavy-ion collision 
processes. 
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APPENDIX A: THERMAL FIELD THEORY 



In this Appendix, we describe the CTP formalism in flat spacetime in order to compute the effective action for a 
scalar field in thermal equilibrium p5| , ^3| -[55|] . Using the CTP formalism one can deal with non-equilibrium situations 
because we do not have to assume anything about the final state of the system in contrast to other formalisms of 
quantum field theory at finite temperature. [ p6| , p7p9| ]. We begin with a brief general description of the CTP formalism 
for an arbitrary initial state. Then, we explain with more detail the particular features of the formalism for an initial 
state in thermal equilibrium. 



1. Arbitrary initial state 

The CTP formalism is a functional technique to obtain expectation values of products of quantum field operators 
and dynamical equations for the mean field which are real and causal |47|j4^ . For simplicity, consider a scalar 
field <I>(x) described initially by the time independent Hamiltonian H[^{x)]. Suppose that an external classical source 
J{x) is turned on and the perturbed Hamiltonian assumes the form H + J{x)^{x). If the initial quantum state of the 
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field is defined by the density matrix operator p, we can define the CTP generating functional, Zp[J], as the following 
statistical average over the initial quantum state [p^ , ^5|j3^ , ^ , 

Zp[J] = Tr{pTce'/c'*"^ '^(^)^(^)}. (Al) 

In the above definition (pix) is the scalar field operator which evolves in time under the unperturbed Hamiltonian H, 
n is the number of spacetime dimensions, Tc is the time ordering operator defined over a contour C in complex time 
which goes forward in time and then back to the initial point (see Fig. 1) and the integral symbol means complex 
integration along this path. Note that in this picture the time evolution of the initial state is determined by the 
external source J{x) alone. Our main concern is to apply this formalism to a system initially in thermal equilibrium. 



Im{t) 




FIG. 1. Integration contour C in the complex time plane for the CTP formalism of non-equilibrium quantum fields. This 
path can be splitted in two separated segments C+ U C-, the first going forward in time and the second backward. The small 
deviation from the real axis is needed to ensure a well defined path integral representation of the generating functional. 



By differentiating with respect to the external source and then taking J — this functional generates the expectation 
value over the initial quantum state specified by the density matrix p of time ordered products of the field operator 
along the contour C. In particular, the two-point function is given by. 



1 



S'Zp[J] 



Zp[0] i5J{x)i5J{x') 



j=o 



Tr{p Tc^{x)^ix')} 
Tr{p} 



(A2) 



To write an integral representation of the generating functional we have to introduce a complete basis of eigenstates 
{|(/9,ti)} of the field operator 4>{x) at the initial time U, i.e. (j){ti,x)\Lp,ti) — Lp{x)\ip^ti) where ^{x) is the eigenvalue. 
This representation for Zp [J] is 



(A3) 

The last term in the integral may be interpreted as a transition amplitude from the initial state {(p^ti) to another initial 
state \ip',ti) following the close time path C, and then, it can be expressed as a path integral with the appropriate 
boundary conditions. 



Zp[J] = / %]d[vp'] {ip,U\p\^',U 



J (C[<k(x)] + .J(x)4,(x)) 



(A4) 



In the above equation, C[(j){x)\ is the Lagrangian density corresponding to the unperturbed Hamiltonian of the scalar 
field theory. Now, as usual, we introduce the generating functional Wp[J] = —i \nZp[J] that generates the connected 
part of the n-point functions and define the CTP effective action as its Legendre transform. 
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r^^^[0] ^ Wp[J]- d-xJix^ix), (A5) 



c 



where (f>{x) is the expectation value of the field (f>{x) over the initial state in the presence of the external source J{x), 

Note that (f) becomes the mean value of the field if we substitute J = 0. The CTP effective action is the generating 
functional of one-particle irreducible graphs and contains all the quantum corrections to the classical action. Since we 
have only needed information about the initial state of the system a real dynamical equation with causal boundary 
conditions can be derived from this effective action . 



Im{t) 




FIG. 2. Integration contour in the complex time plane for the CTP path integral useful for a system with an initial state in 
thermal equilibrium. 



2. Initial state in thermal equilibrium 

For practical purposes, we are going to describe in some detail how to obtain the Green functions and the effective 



action for a quantum system with an initial state in thermal equilibrium at temperature /3 ^ [£4 38 55[| . Consider the 
Lagrangian density for a free massive scalar field 

mx)] = -i(a^(/.a''0 + m2</)2). (A7) 

The initial thermal state of this system is described by a normalized density matrix of the form — e^^^ /Tr{e~^^}, 



where H is the Hamiltonian operator corresponding to (A7) at some initial time U; i.e., the unperturbed Hamiltonian 



of the system before the external source J{x) is connected. In this case, the path integral representation of the 



generating functional given in Eq. (A4) can be written as 



1 r r 



Tr{e-'3^} 



In deriving this expression, we have used the fact that the density matrix operator pfs can be seen as a time transla- 
tion operator in the complex plane because the unperturbed Hamiltonian is the generator of time evolution for the 
eigenstates of the field operator (j){x). In other words, since the field transforms as (jiitjx) = e'*^0(O, i:)e~'*^ and the 
eigenstates evolve in time simply as \(p,t) — e^*^\tp,Q), the density matrix generates a translation in complex time 
t t — i[3. The above representation for Zp[J] may be written in a more compact form as 
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The path integral is defined over all possible field configurations along the time path C with initial boundary condition 
(j){ti) — (j){ti — ij3) and the normalization term M includes all the factors which are independent of the external source 
J{x). We would be also able to write the above compact expression for the path integral along the time path in 
the complex plane of Fig. 2. But, the contribution to the path integral along the segment C/5 has no dynamical 
consequences in our case because the external source is zero along this vertical segment. In the usual path integral 
approach to thermal field theory [^ , ^7|j4g[| (see also [Q), the external source is also non zero along Cp by definition 
in contrast with our approach |38 4^ , |55|| . However, it can be justified that both are consistent asymtotically; that is. 



when ti —00 and — » 00. We are going to assume this condition from now on. 

To perform the path integral of the generating functional it is convenient to make the change of field variables 
(f){x) 4'{x) ~ J^d"x' G'^{x — x')J{x'), where the propagator G'^{x — x') satisfies the differential equation, 

{Uc-m^)Gfi{x~- x') = 5c{x-x'). (AlO) 

The subscript C in the d'Alambertian Dp and the delta function 5c means that these objects are defined along the 
complex path C . The subscript /3 in is written to remind that this propagator has to be constructed with the 
appropriate thermal boundary conditions. If we redefine the normalization factor J\f by 2^/3 [0] in order to include a 
new path integral independent of J(x), the change of variable yields 

Zp[J] ^ ^^[0] g-i/o''"-'^"-' (All) 



It is clear from the thermal version of (Al) and the above expression for the generating functional that the propagator 
G'^{x — x') can be interpreted as the thermal average of two time-ordered operators along C 

G'fj{x-x') = ~i{Tc<l>{x)cl>{x'))p. (A12) 

Above we have introduced the shorthand notation {A)fj = Trlp/^A} for the trace of a general operator A with respect 
to the thermal density operator. It does not take more additional work at this point to show that the Wick theorem 
also applies at finite temperature fs^ , ^ . Let us find an explicit expression for the thermal propagator. Using the 
step function 9c ^ defined also along the contour G, one can decompose the propagator as 

G^(r-r';x-f') = Ocir - t') G+{t - t';x - x') + Ocir' ~ t) G^{t - t';x - x'), (A13) 

where we have denoted the time variable along the path G by t and defined the propagators G^(r — t' : x ~ x') = 
G~^{t' — t'x' — x) = —i{4>{T]x)4>{T']x'))f}. The trace nature of the generating functional and the Hamiltonian time 
evolution of the field operator 0(x) impose some relations on certain propagators. For example, they lead to equalities 
like G^ (r; x) — G^ {r — if3; x) which are generally known as Kubo-Martin-Schwinger (KMS) relations ||l^,|ll| . In fact, 
these conditions are a direct consequences of thermal equilibrium and are necessary to preserve the periodic boundary 
condition in the change of field variable previously performed in order to simplify the path integral representation 
of the generating functional [ p8|j59| ]. The solution of ( AlO| ), compatible with the KMS relation, may be expressed in 
terms of the thermal function nfj{p°) — {e^P° — 1)^^ as 

/fin ^ 
_JL e'f(-- ) [6c{t - t') + np{p°)] {-2m) sg{p°) S{p' + m'). (A14) 

Now, for simplicity in actual computations, it is convenient to introduce a matrix notation for the thermal propagator 
Gp which takes into account the segment of the path, G+ or G_, where the time variable is choosen to lie. We have 
four different possibilities: 
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• if T, r' S C+ or r, r' S we may write, respectively, 

Tl 



, .--27rin0(|p°|)(5(p^ + m^) 



(A15) 



if r e C+ and r' G C_ , or vice versa, then respectively, 



Since the time variable in the above equations takes values on the real axis alone, it is very easy to see that all these 
propagator components may be expressed as thermal averages 

gI+{x-x') - -i{T4>{x)4>{x'))fi 
GtA^-x') ^ -i{f4){x)4){x'))p 

gP_+{x-x') = -i{cj,{x)cj,{x'))p = g1_{x'-x), (A17) 

where T and T are the usual time and anti-time ordering operators, respectively. If we define in the same way the 
external sources J±{t, x) = J{t, x) for r e C±, we can write the generating functional in matrix form 

In the exponents, the time integrals range from —oo to +oo, the subindixes a, b take values +, — and the transposed 
source vector is defined by Ja{x) = ( J+{x),—J^{x) ). The minus sign in the second component of the source 
vector is a relic to the fact that originally the integral for the time variable along the segment C_ was from +oo 
to — oo. If we were interested to specify initial vacuum boundary conditions (po = |0,m)(0,m|) instead of thermal 



equilibrium conditions we could still use (|A18|) substituting the thermal propagator G^^ by the vacuum propagator 
G°fc = hm/3_oo eft (see, for example |37|j4^,|55t ) . 

Finally, we are interested to obtain the CTP effective action at finite temperature in terms of the propagator. 
Following [Q, where the derivation of the CTP effective action in the case of an initial pure vacuum state was 
explained in some detail, but taking into account the particular boundary conditions for the thermal case, we arrive 
to the expression (see also ]55| , p3| ) 

T^CTp[4>a] = 5[0a]-^Tr{lnGfj, (A19) 

where S[(t)a] is the classical action for the scalar field evaluated at the expectation value of (t>a{x) in the presence of 
the external source Ja{x). 

3. CTP effective action for two fields in the semiclassical approximation 

We now derive the CTP effective action at finite temperature for two fields in the semiclassical approximation in 
order to generalize Eq. ( A19| ) and make it of direct applicability to our problem. For simplicity, consider two different 



scalar fields described by the general Lagrangian density 

C[ci,{x)]^C[iP{x),<j,{x)], (A20) 

where the second term is assumed to be quadratic in the field il^ix). In our case, the first part, C[(j){x)\^ will represent 
the Lagrangian density of the gravitational field, and the second part, the Lagrangian density of a thermal scalar field 
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coupled to gravity. Since, we are only interested in the dynamics of the gravitational field we perturb the system by 
an external source coupled only to (t>{x). The CTP generating functional for this system will be 

^iW(i[J±] ^ qy^^-^jy^^^ J^d"x (C[<p{x)]+C[-ip(x),4,{x)] + ,J(x)4,{x)) (A21) 



where C is the time integration path of Fig. 1 and the field ^{x) has the appropriate thermal boundary conditions. 
If one of the fields is assumed to evolve classically, in our case (j>{x), the above CTP generating functional may be 
approximated by 

giWi3[J±] ^ Ic'^"^ iC[cj>{x)] + J{x)(j)(x)) / pj^j J^d"x C[4>(x),<p(x)] (A22) 



where (t>{x) is the classical field written in terms of J{x). The remaining path integral is in fact a functional of 4>{x) 
alone and can be identified with the influence functional of Feynman and Vernon |ll5|| . Since C[il!{x), (l>{x)] has been 
assumed to be quadratic in ip{x), the path integral is Gaussian and may be performed exactly. In matrix notation we 
have 

exp(-lTr{lnGfj0±]}), (A23) 

where G^f^[(j)±] is the corresponding matrix propagator with thermal boundary conditions of the Lagrangian density 
C[ip{x),(j){x)]. When this result is substituted in the CTP generating functional and the Legendre transform with 
respect J{x) is done, we can finally write the thermal CTP effective action as 



C 

This generalized expression for the CTP effective action justifies the use of Eq. (^. 



r^Tpfe] = S[4>+] - S[4>-] - -rr{lnGfj0±]}. (A24) 



APPENDIX B: TENSOR BASIS 

The complicated tensor structure of the thermal CTP effective action may be simplified if the following tensor basis 
is used 

Tp"^{X,Y) = X^X^^Xf^X" 

Tf^"''(A,y) = r|^"'r|"'^ 

T^'''"''(A,y) = A^A'^77"'' + A"A V" 

T^''^"''(A,y) = A''(y''?7"^ + r^77"'^) + A"(r''r;'3^ + y^ry'5'') 

+y'3(A'^77"^ + A'^Ty"'') + y"(A''/^ + A'^Ty'^'^) 
T^'''"''(A,y) = r'^A^A'^A^ + r^A^A'^A'^ + r'^X^A^A'^ + r^A'^A'^A" 

T^'''"''(A,y) = yA'y^x^A'' + y^y^A^A'^ 

T5'o'''"''(A,y) = (y'3A° + y"A'')(y''A^ + y'^A'') 

rp/ji/.Q/s^j^ y-j _ Y^Y^^Y'^ X" + Y"Y^Y'^ X^ + Y"Y^Y'^ X^ -\-Y'^Y^Y^X'^ 

T^2'"''(^.^) = y^y^y^y^ 

T5'3^^"''(A,y) = yA'y^^a/S _^ y"y/3^M^ 

T5'4''^"''(A,y) = (y''A^ + y'^A'')7y"'' + (y''A" + y"A'^)7y^''. (Bl) 
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This basis of 14 tensors with four indixes is constructed with two different four- vectors X^^ and Y'^. Note that in 
order to maintain the symmetries of the effective action each member of the basis has to be symmetric under the 
interchanges ^ <-> i^, a ^ /3 and [^xv) ^ {ctfi) 0]. In terms of this tensor basis, the tensor T^'^'°'^{q^ fc), introduced in 



Sec. [IB, is defined by 



i:^''''^P{q, k) ^ Tf ■"'^(g, fc) + . 9) + q^Y + ik\k ■ q) + q^] + C^fc^j T^'',"/5(g^ fc) 

\[{k ■ q) + q^] + ik'^ Tf '"'^(g, k) + iT^'"''!?, k) + ^T'^'^-'^^q, k) 
+ k) + l^.-^tr'iq, k) + e-^tr'i^^ k) 

{\[{k-q)+ q'] + TZ-'^^q, k) - \l\[{k ■ q) ^ q^] + ^k^X T^.^^^^q, k). 



(B2) 



APPENDIX C: SOME COEFFICIENTS AT HIGH TEMPERATURE 



The 14 coefhcients Ili(r) needed to write the leading contribution at high temperature of the kernel h'^ ' ^ {k) in 
terms of the tensor basis of Appendix B are 



Hi 

H2 

H3 
H4 
H5 
He 

H7 

Hg 
Hg 
Hio 
Hii 

H12 
H13 

Hl4 



1 1 



1 



6 24 



3 12 



8 
5 

24' 



8 



3 12 24 8 

24 8 

--if2_1^4_ 5^6^ 

12 24 8 



12 24 
1 7 



35 



35 



3 12 



24 



1 



= + —K^ - -K^L + -K^L 

12 24 2 8 

= _i_2^2 35^4 15^4^ 35^6^ 
6 3 24 4 8 



1 2 



35 



15 



35 



= + —K^ K^L + —K^L 



6 3 
1 35 
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4 24 2 8 y 

12 24 2 8 

----if2_ 5 4 
12 24 8 



(CI) 



where K^^ = k^^ /\k\ = (r. A;) and the function L is given by 

L(r) 



^Inf^ 
2 \r-l 



(C2) 
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